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Abstract 

Theoretical studies of jet stopping in strongly-coupled QCD-like plasmas have used gauge-gravity 
duality to find that the maximum stopping distance scales like i?^/^ for large jet energies E. In 
recent work studying jets that are created by finite-size sources in the gauge theory, we found an 
additional scale: the typical (as opposed to maximum) jet stopping distance scales like (EL)^/^, 
where L is the size of the space-time region where the jet is created. In this paper, we show that 
the results of our previous, somewhat involved computation in the gravity dual, and the {EL)^/^ 
scale in particular, can be very easily reproduced and understood in terms of the distance that 
high-energy particles travel in AdSs-Schwarzschild space before falling into the black brane. We 
also investigate how stopping distances depend on the conformal dimension of the source operator 
used to create the jet. 



I. INTRODUCTION AND RESULTS 



Various authors p}M] have made use of gauge-gravity duahty to study the stopping dis- 
tance of massless, high-energy jets in a strongly-coupled plasma of A/'=4 supersymmetric 
Yang Mills theory (with and without the addition of fundamental-charge matter). All have 
found that the furthest that such a jet penetrates through the plasma scales with energy as 
E^^'^. Most of these methods specified the initial conditions of the problem in the gravity 
description of the problem, and it is not completely clear exactly what these initial con- 
ditions correspond to in the gauge theory. However, we recently showed j3] one possible 
way to set up the entire problem directly in the gauge theory, only then translating to the 
gravity description using the conventional elements of the AdS/CFT dictionary. We specif- 
ically studied jets that carried R charge, and we measured how far that charge traveled 
before stopping and thermalizing. Though we did find that the furthest charge would travel 
through the plasma scaled as E^^^, we also found that, on average, almost all of our jet's 
charge stopped and thermalized at a shorter distance that scales as (ELY^'^, where L is the 
size of the space-time region where our jet was created. Fig. [T] shows a qualitative picture 
of our result for, on average, how much of our jet's charge was deposited as a function of 
distance traveled x^. (Our convention here is to write 4-dimensional space-time position as 
x'^ and take our jets to be created near the origin, traveling in the x^ direction.) Between 
the (ELY^^ scale and the E^^^ scale, the distribution falls algebraically like {x^)~^ for jets 
created by the source used in Ref. We will work in units where 27rT = 1, and in those 
units the specific formula we derived for Fig. [T] was 

P^oHx^) - 2^^^ *(-^) for x3 « i?V3 (i.ia) 



and 



Prob(x') ~ 4^^^ ^(0) exp f-^i^l for x^ > E^/\ (1.1b) 

E \ E^'-^ J 

where "^{y) is a source- dependent function that suppresses \y\ ^ 1, causing suppression of 
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X" <ti (ELY^^ above. The c's are constants given by 



-(^. (1.2) 

ci ~ 0.927, and Ca ~ 3.2. 

The calculation that produced (11. ip was long and not particularly enlightening as to the 
origin of the {ELY^^ scale. The purpose of the current paper is to show how that scale, 
and then the precise result (II. lap for the case x^ <^ E^^^, can be derived from a very simple 
calculation of how far a classical massless particle travels in AdSs-Schwarzschild space before 
falling into the black brane. In the process, we will learn more about exactly what feature 
of the source determines the {ELY^'^ scale. We will see that it is not directly the size 
but the typical "virtuality" of the source that matters (where = q^rj'^'^qu is squared 
4- momentum) . 

Our analysis of the distance traveled by falling particles in AdSs-Schwarzschild will be 
essentially the same as an earher analysis by Gubser et al. [l[ and Chesler et al. [3|, who 
used it in a discussion of the falling endpoint of a classical string. The difference here will 
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FIG. 1: The average deposition of charge as a function of for jets created by the source described 
in Sec. HVAland in Ref. fl. 



be one of context and application: Our analysis of jets [4i] does not involve classical strings, 
and we will use the falling particles to explain the {ELY^^ scale. 

In our earlier work we created the jet by turning on a small-amplitude source whose 
space-time dependence had the form 

source(x) ~ e*'' ''AL(x) (1.3) 

of (i) a high-energy plane wave e^^'^ times (ii) a slowly varying envelope function Al{x) that 
localizes the source to within a distance L of the origin in both space and time. We took k 
to be light-like: 

F = (E,0,0,E). (1.4) 

In addition, for the sake of simplicity, we took the source to be translation invariant in the 
two transverse directions. So, for example, 

A^(a;) = e-^(^We-^(^W. (1.5) 

The Fourier transform of the source (11. 3p is non-negligible in the region of momentum space 
depicted in Fig. [2^: a region centered on k with width of order L~^. We take ^ E. 
Note that this source covers a range of values of g^, from to order ±E/L, and the typical 
size of is order E/L. 

In the gravity description, this source causes a localized perturbation on the boundary of 
AdSs-Schwarzschild space-time, which then propagates as a wave into the fifth dimension, 
eventually falling into the black brane horizon. The analysis of jet stopping in Refs. [2I, 0] 
was based on the analysis of such 5-dimensional waves. 

Now imagine instead a source where k is slightly off the light-cone, 

F = (E + e,0,0,E-e) (1.6) 

with e <^ E, and where the envelope size L is wide enough that the picture in momentum 
space looks like Fig. [2]d instead of Fig. |2^, with the spread 1/L in momenta small compared 
to e. In this case, the of the source is approximately well defined, with ~ A;^ ~ 
—AEe. We will show that in this case the wave created by the boundary perturbation 
is localized into a small wave packet, whose motion may be approximated by that of a 
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FIG. 2: Qualitative picture of momenta contributing to the source (|1.3p used to generate jets (a) 
for the calculation originally used to find with <^ E, and (b) in the case <^ e <^ E of 

(jl.6p . Figure (c) depicts (a) as a superposition of distributions of type (b). The cells in (c) that are 
extremely close to the light cone cannot be treated in particle approximation, but the contribution 
of cells that can be treated so dominates when L <C the maximum stopping distance scale E^/^. 



classical, massless particle which starts at the boundary, traveling in the direction, with 
4-momentum proportional to q. The trajectory of such a particle is shown qualitatively in 
Fig. [Sti. By a simple calculation, we will find that the particle falls into the horizon after 
covering a distance 



r\ , 

"stop — I _ 2 



I 12 \ 1/4 / 1/4 



V2\-qy 2 



(1.7) 



where the constant c is given by (11.21) . As measured by boundary time x^, the particle takes 
an infinite amount of time to fall into the horizon. As it gets closer and closer to the horizon, 
the boundary distortion that the particle creates (see Fig. [Sjo) becomes weaker and more 
spread out, which corresponds to charge diffusion in the boundary theory after the jet stops 
and thermalizes. This qualitative picture is similar to the qualitative picture of the effects 
of a classical string falling into the horizon given in Refs. [H, |6|.^ 

Note that the stopping distance (II. 7p only makes sense for < (i.e. e > 0). The q^>0 
components of a source do not create an excitation of the system that persists after the 
source turns off and so are not relevant [2I, 0] . 

Now consider the original source of Fig. as a superposition of sources like Fig. IJJd, as 
depicted in Fig. Since sources with different values of e have different stopping distances 
(II. 7p . we might guess that the different pieces of this superposition do not interfere and 
so the source of Fig. simply produces a distribution of stopping distances, weighted by 
independent probabilities that the source produces a jet with a particular g^. That is, 

Prob(x^) ^ j %2)P(g2)(5(x^-xLp(g^)), (1.8) 

where V{q^) is the probability density for the source to produce a jet with a given g^, and 



See in particular the discussion surrounding Fig. 2 of Ref. j6|, which inspired our Fig. |31d. See also Ref. 
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(a) (b) 



FIG. 3: (a) A classical particle in the AdSs-Schwarzschild space-time, moving in the direction 
as it falls from the boundary to the black brane in the fifth dimension u. (b) The presence of the 
particle (the large dot) perturbs the boundary theory in a manner that spreads out diffusively as 
the particle approaches the horizon for ^ oo. 



where the stopping distance xl^^^^q'^) is given by fll.7p . We will verify that this formula 

precisely reproduces the -C E^^^ case (11. lap of our previous result. 

We can now see where the {ELY^^ scale comes from. It is only time-like source momenta 
< that produce jets. The typical value of time-like for the source of Fig. [2^ is 

g2 ^ —E/L, corresponding to e ~ L~^. Putting this into (11. 7p . the typical stopping distance 

in this case is therefore 

<pical ~ (^^)'/^ (1-9) 

Note that it is the of the source that determines the stopping distance, and that the 
typical value of is determined by L in the case of Fig. |2^. 

The estimate (II. 9p of the stopping distance ceases to make sense if the size L of the source 
becomes as large as the stopping distance itself. This happens when 

L-xl^^^{ELy/\ (1.10) 

which gives 

4op~^'/'- (1-11) 

We will see later that this is precisely the case where the wave packet in AdSs-Schwarzschild 
can no longer be approximated as a particle. The moral is that the simple particle picture 
gives us not only the {ELY^^ scale but also, simply by estimating where it breaks down, the 
E^^^ scale as well. 

In the next section, we will briefly review the trajectories of massless particles in AdSs- 
Schwarzschild and derive the corresponding stopping distance (II. 7p . In section ilTTI we discuss 
the conditions for being able to approximate the 5-dimensional wave problem with particle 
trajectories and verify that they apply in the case of interest. Then we use the particle 
picture in section |IV] to simply reproduce our original result (II. lap for charge deposition for 
x^ ^ E^^^. In section lYl we generahze our results to jets created by other types of source 
operators than those originally considered in Ref. 0]. We will see that Fig. [T]is modified to 
Fig. m Finally, we offer our conclusions in section IVII 
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FIG. 4: The probability distribution of jet stopping distances for scalar or transverse BPS sources 
with conformal dimension A. Fig. [T] corresponds to A=3. The scales a^typicai ~ {EL/\/~KY/^ and 
^max ~ (-B/A)^/^ indicated along the vertical axis assume that A is held fixed when taking the 
limit of large energy E (as well as large coupling g^Nc and large Nc). The parametric scaling with 
A indicated for x^^^^^^.^ assumes a Gaussian source envelope (jl.Sp . but the other features shown in 
the figure are independent of the details of the source envelope. 



As an aside, some readers may be curious how the stopping distance scales {ELY'^ and 
^1/3 generalize to other dimensions. On the gravity side, it is very easy to generalize the 
results of this paper to different space-time dimensions d of the boundary, but it is less 
certain what strongly coupled field theories these classical gravity theories correspond to. 
(See Refs. [3] for proposals.) Ignoring the question of interpretation, we show in Appendix 
|A]that {Eiy/^ and E^^^ generalize to (EL)^'^-^)/^'^ and E('^-2)/(rf+2) respectively for d>2. 

In this paper, we will use the convention that Greek indices run over the 4 space-time 
dimensions [fi = 0,1,2,3) of the boundary theory and capital roman indices run over all 
five dimensions (/ = 0, 1, 2, 3, 5) of the AdSs-Schwarzschild space-time. The symbol will 
refer to the squared 4-momentum, = q^rj^^qu = — + When we use light-cone 

coordinates, our conventions will be 

V^ = V^±V\ V± = \V^ = \{y^TV^) (1.12) 

for any 4-vector V . Throughout this paper, the adjective "transverse" will refer to the 
spatial directions 1 and 2 orthogonal to q.^ 



II. REVIEW OF FALLING MASSLESS PARTICLES 

Null geodesies in a 5-dimensional space with 4-dimensional translation invariance are 
given by (see Appendix [Bll : 

.n^.^) = j^.d.^j^£^,, (2.1) 



^ This is as opposed to the alternative usage of "transverse" to mean all three space-time components 
orthogonal to q^. Note that the spatial directions transverse to q are the same for all q in our problem 
because we take our source (14. 2p to be translationally invariant in the transverse directions. 
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where g is the 5-dimensional metric and qi is a constant of motion for J = 0, 1, 2, 3. We will 
work in coordinates where the metric is^ 



-{-fdt' + dx^) + ^du^ 



(2.2) 



where / = 1 — m^, and R is the AdSs radius. The boundary is at m=0 and the horizon at 
u=l. If we take the 3-momentum q to point in the direction, writing = {—u, 0, 0, \ q\), 
then (12.1 p gives the total distance traveled in falling from the boundary to the horizon 
to be 

du 

(2.3) 



X 



stop 



W' 



where = q^rj^^^qu is the fiat-space square of the 4-momentum. This is the same result as 
Refs. 0,3.^ 

Now let us apply this result to the case \q^\ <^ |qp ~ E'^ relevant to the source of Fig. [2}3. 
For small — g^, the integral of (12. 3p is dominated by small m, and so we may approximate 



■^stop — 



du 



c 

7! 



1^1 



2\ 1/4 



(2.4) 



which gives (II. 7p . 



III. WAVE PACKETS AND GEOMETRIC OPTICS 

In this section, we will discuss the conditions necessary for making the particle approx- 
imation. A wave packet behaves like a particle when it is wide enough to contain many 
phase oscillations of the field yet small enough that the properties of the background do not 
vary significantly across its width, as depicted in Fig. |5]d at a particular moment in time. 
We can arrange such a width provided the background properties do not vary significantly 
over one wavelength of the phase oscillation. This is the geometric optics limit, which we 
referred to in our earlier work [J] as a WKB approximation. To check the geometric optics 
limit, one may focus as in Fig. [5}:; on a wave with a single, generic value of q typical of the 
wave packet, and investigate how much things change over one phase oscillation. 

To assess whether a wave packet is adequately particle-hke to use a particle-based cal- 
culation of the stopping distance, it will be helpful to understand what the important scale 
for u is in determining the stopping distance (12. 4p . The integral in (12. 4p is dominated by u 
of order 




Our formulas in this paper would be a little tidier (fewer square roots) if we worked with the coordinate 
z = instead of u. We will stick with u in order to facilitate coirmarison with our previous work Q]. 
^ Our (|2.4p corresponds to the first part of Eq. (5.3) of Gubser et al. P|, where their y is our -y/u, their zr 
is 2 in our units 27rr=l, their pi/po is replaced by the qs/qo of the momentum Ou typical of our source, 
and their yuv is set to zero. It also corresponds to Eq. (4.28) of Chesler et al. [3|, where their u is our 
2^/u, their Uh = 2 in our units, their ^ is replaced by our qo/qi, and their is set to zero. 
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A "l/qs A ~"l/,5 A ~"l/<,j 



(a) (b) (c) 

FIG. 5: (a,b) A snapshot in time of waves in the fifth dimension u for times after the boundary 
source has turned off but early enough that n < n^, ^ 1 (that is, before the wave gets very close to 
the horizon), (a) shows the type of wave generated by a localized source that superposes a range 
of values such as Fig. [2^. (b) shows the wave packet generated by a source with approximately 
well-defined such as Fig. [2)o. (c) shows a single 4-momentum component, corresponding to a 
single, definite value of 4-momentum q^. Case (a) matches Fig. 7a of our earlier paper [J]. 



The relevant question is then whether the background varies significantly across one phase 
oscillation for u ~ u^. By fll.7p . the distances xl^^^ < E^/'^ relevant to jet stopping correspond 
to e > E~^/^ in the particle picture and so to u^, < E~'^^^ <^ 1. So we may focus on small u 
in what follows. 



A. Geometric Optics 

For a massless 5-dimensional field with definite 4-momentum q^, the exponential in the 
WKB approximation is 




where the are constant and q5{x^) is determined by the 5-dimensional massless condition 

qig^'^qj = o, giving^ 

(l5{x^) = 955{-qf, 9'"' Qu)- (3.3) 

For the metric (12. 2p . this is 

. i^^E. (3.4) 

For the geometric optics limit, we need the wavelength to vary insignificantly over one 
wavelength.^ For the important values m ~ ^ 1 of m, the wavelength X{u) ~ l/qs{u) in 



^ The null geodesies (|2.ip can be expressed in terms of 55 as x'^{x^) ~ — J dx^ {dq^ / dqf^) . This particle 
formula is simply the saddle point condition dq^ [iq^x'^ + * / dx^ qb{x^)] = with respect to for the wave 
(|3.2p . Also, / dx^ qsix^) was referred to as the WKB exponent 5 in Ref. Q], where various expansions of 
the integral may be found. 

^ This condition can be phrased in an z^-reparametrization invariant way as V5(l/(75) ^ 1, which is 
(955)^^ dsKgssy / qs] <C 1. We give more detail in Appendix [Cl in the specific context of the particular 
type of source operator that we used in our original calculation. 
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the fifth dimension satisfies this condition if 



u^q^iu-,) > 1, (3.5) 

which, using fl3.4p . is^ 

> (3.6) 

-r 

Using the size (13.11) of m^, this condition is {—(fY I\q.\'^ ^ 1, or equivalently 

e > E-^^'^. (3.7) 

Referring to the stopping distance (11. 7p . we then see that the approximation of geometric 
optics, necessary for a particle interpretation, breaks down unless 



X. 



stop « E''\ (3.8) 



Within the context of our approach to jet stopping, a proper analysis of what happens at 
distances > E^^^ requires a wave rather than particle description of the problem, as in Refs. 
[21, 0]- The wave analysis gives exponential fall-off for propagation beyond i?^/'^, as described 
by (HlbD. 

More on the geometric optics approximation can be found in Appendix O 



B. Wave Packets 

The geometric optics limit for (13. 8p allows us to create localized wave packets. Here we 
will see how wide those wave packets are in u for sources of the form of Fig. [2]d. We are 
primarily interested in the case where the center of the wave packet is at the critical scale 
M ~ in the fifth dimension. However, the presentation will be a little more straightforward 
if we first make parametric estimates for earlier times, when the center of the wave packet 
is at M <C M^, and then extrapolate those parametric estimates to m ~ m*. 

For u -C M*, the distance traveled (12.11) for a massless particle is 

x\u)= r^^=^2,l^^^l-. (3.9) 
Turning this around, the location of the particle in the fifth dimension is 

(3.10) 

When the particle is replaced by a wave packet, there are two sources of uncertainty. The 
size L of the source introduces an uncertainty in the initial position of the excitation of 
Ax^ ~ L. It also introduces an uncertainty in the 4-momentum q, and so in e, of Ae ~ 1/L 
as in Fig. From (I3.10p . the combined uncertainty Am in u is then of order 




In the language of Ref. |J|, the condition Ui, ^ l/(— q ) is >• Mmatch- 
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Extrapolating this parametric estimate to the case m ~ Mj, of interest, (13. ip and (13. lip 
give 




(3.12) 



The wave packet will then be localized provided (i) L ^ 1/e, as in Fig. |2]d, and (ii) L <^ 
{E/eY^^. By (II. 7p . the last condition is just the condition that L be much less than the 
stopping distance xl^^^. 

IV. REPRODUCING THE DISTRIBUTION OF STOPPING DISTANCES 

Now we will show that the formula (II. lap found in our earlier work [4] for the average 
distribution of charge deposition can be understood as a convolution (II. 8p of the particle 
stopping distance with the probability density V{q^) for the source to create a jet with a 
given q^. 

We will later give in section |V]a very general argument, based on dimensional analysis, for 
determining the 'P{q^) associated with different choices of source operator. This argument 
will also require a discussion of massive fields and massive particles in the gravity dual. For 
the moment, we will be less general, stick to the specific type of source operator that we 
used in previous work, and show how to extract 'P{q^) from a result for the average total 
charge produced by the operator. Readers who would prefer to just see the more general 
argument may skip section IIV Al below and instead wait for section IV Bl 



A. Extracting T'{q^) from results in the literature 

In Ref. we used a source involving R-current operators j^. Specifically, we modified 
the 4-dimensional gauge theory Lagrangian by 



with a localized background field 



(4.1) 



Ak-x 



L 2 



+ h.c. 



Al(x), 



(4.2) 



where Ma is an arbitrarily small source amplitude, e is a transverse linear polarization, and 
r* are Pauli matrices for any SU(2) subgroup of the SU(4) R-symmetry. We then measured 
the response {j^^^^{x)) of the R charge current associated with r^/2. The gravity dual to 
the R charge current operators is a massless 5-dimensional SU(4) gauge field. We chose our 
source (14. 2 p to be translationally invariant in the transverse spatial directions {x^,x^) to 
simplify the calculation. In what follows, we will refer to the R charge associated with r^/2 
as simply "the charge." 

For an arbitrarily small source amplitude A/a, the source will usually have no effect at all 
on the system. On rare occasions, with probability proportional to A/"^, the source creates 
an excitation (in our case a "jet") with the same quantum numbers as the source operator. 
In our case (14.21) . that means it creates a jet with total charge equal to 1. The creation 
of an excitation with different quantum numbers would be even higher-order in A/^ and so 
negligible. Since an excitation (if any is created) has charge 1, the average total charge Q 
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created by the source is then equal to the probabihty of the source creating a jet. We will see 
that from the previously calculated result for Q we can then extract the probability density 

In Ref. 1^,^ we showed how to use the field theory Ward identity to make a simple 
calculation of the average charge Q. Here we will just quote the result, which was 



Ik J W) 



I -j^J{-q')\q'\Mq-k)\\ (4.3) 



where g^Q = Att/Nc {Nc^oo is the number of colors) and 6 is the step function. The 
-^L{q — k) in 04.31) is simply the Fourier transform of the x-dependence fll.3p of the source. 
Our source only has support for q = k + Aq with Aq small compared to k, in which case 
~ 4:Eq^. Eq. f l4.3p may then be approximated as 

^ SnEAfl C2dq+dAq_^. ,, ,,7(2), . s|2 

Q- / ^l^_^0(-g+)|g+||A^^(g+,Ag_)| , (4.4) 

~ ('2') 

where A.i is the two-dimensional Fourier transform 



of the source envelope and V± is the area of transverse space {x^ and x^). Since (14. 4p involves 
an integral over g+, and since Q is the total probability of creating a jet, it is natural to 
interpret (14. 4p as giving a probability density 

167lE^^l^^ f dAq_ „ iIa(2)/ a„ m2 



P+(g+)^ / —i-e{-q^)\q^\\A'^>{q^,Aq^)\' (4.6) 

fi'SG J l^^J 

for producing a jet with a given value of g+. If a jet is produced, the probability distribution 
for its g+ is then the relative probability 

Here we've put a bar over the argument of P+ just to distinguish it from the g+ integration 
variable in the denominator on the right-hand side. Now package the source dependence 
into the definition 

AL^ Jdq+ dq^ 0{-q+) | ^ (g+, g_) | 





1?+ 






\f (g+,Ag_)| 


2 


J dq+ dAq^ 6 






\\Kf{q+,Aq.)\' 



This is the definition we made in our earlier work for the \l/ that appears in the charge 
deposition result (II. ip . With this definition, 

P+(g+) = AL^ ei-q+) \q+\ ^(q+L). (4.9) 



^ See specifically Appendix A of Ref. Q| . 
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B. Using P(q^) to get Prob(a;^) 



Since ~ 4£'g+, the probability distribution V{q^) for is related to the probability 
distribution V+{q+) of (14.91) for by 

However, at this point it will be easier to just stick with g+ and V^. 

Note that in each cell of Fig. the typical value of g+ is just what we have previously 
called — e for that cell. The distribution (11. 8p of stopping distances based upon the picture 
of massless falling particles in 5 dimensions can then be written as 



Prob(x^)~ J rfeP+(-e)5(x3-xLp(e)), 



(2x3)9 V (2x3)4 

As promised, in the case x^ <^ E^^^ where we have argued that the particle picture should 
work, this exactly reproduces our earlier result (II. lap that came from a full, much more 
complicated calculation. 



V. MASSIVE PARTICLES IN 5 DIMENSIONS 

In the preceding sections, we have assumed that the 5-dimensional bulk field dual to 
the source which creates the jet is massless, such as the 5-dimensional gauge field dual to 
R current operators. One may wonder what results change if we choose different types of 
source operators that are instead dual to massive bulk fields. In this section, we will see 
that the basic qualitative picture of Fig. [T] of the distribution of stopping distances remains 
the same, except that the exponent of the (x^)"^ power- law tail changes, depending on the 
conformal dimension of the source operator. 

When we wish to make contact with a particular example, we will for simplicity restrict 
attention to scalar BPS (e.g. chiral primary jsf) operators.^ In that case, the mass m of the 
5-dimensional field is related to the scaling dimension A of the operator by [9| 

{RmY = A{A-d), (5.1) 

where R is the radius of AdSs, and d=A is the dimension of ordinary space-time. The 
possible values of A are bounded below by 

A>^. (5.2) 



Examples of scalar BPS operators include the Lagrangian density and the symmetrized trace 
tr((/)(-jj • • • 0^^)), where (/)2, and 03 are the three complex scalar fields of A/'=4 supersymmetric 
Yang Mills. 
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In what follows, we will hold m fixed when we consider the limit of large jet energy E. 
However, it will be interesting to consider the case A 3> 1 (i.e. Rm 3> 1) in addition to the 
case where A is of order one. We will find that the typical and maximum stopping distances 
decrease for larger A. 

We will study the propagation of excitations of a massive bulk field by studying the 
propagation of massive particles in the bulk, similar to the massless case studied earlier in 
this paper. (An alternative discussion directly in terms of a wave analysis is sketched in 
Appendix [Dl) We should emphasize that the term "massive" refers only to the bulk fields 
and corresponding bulk "particles" in our discussion, and so to the conformal dimension of 
the source operators in the boundary theory. We have not introduced any masses in the 
4-dimensional strongly-coupled field theory: the theory is still just A/'=4 supersymmetric 
Yang Mills theory. 



A. Stopping distance of massive particles 

For a particle of mass m in 5 dimensions, the stopping distance integrals (12. ip and (12. 3p 
are modified to 

3 _f / , 5 ^^"g^ _ f (.o-s 

V 1^ 4|(7p J 

As in the massless case, we shall see below that the stopping distance will be dominated by 
M ~ <^ 1. So we will be able to approximate / — 1 above: 

3 f d'^ 

'm(m' - ) - ^ 

For A > c/ (in which case is positive), there is an issue with the lower limit of 
integration in (15. 4p : our classical particle with 4-momentum cannot exist in the region 
where the square root in (15. 4p is imaginary. For u <^ u^. (and focusing on < 0), this 
condition allows for a classical particle when u > Mmin with 

~ (-^^)' 
—4q^ 

How to interpret this? The wave equation is not well described by geometric optics near the 
turning point Umin- However, as long as u^i^ <^ u^, the calculation of the stopping distance 
will be dominated by much larger u (where m is ignorable), and so we may still use the 
particle picture to approximate 

4„„.r , _ .f^L^. (5.6) 



which is the same as the massless particle result (12. 4p . In this respect, the mass can be 
ignored. 
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What happens at m ~ Mmin can be made more concrete by returning to the wave problem 
and looking at the solution to the massive scalar wave equation 

^ -My/^a'^ds^) = (q^g'^q. + m')^ (5.7) 



in the limit m <^ -u^ <^ 1. In this limit, one is close enough to the boundary that AdSs- 
Schwarzschild is approximately just AdSs, and the equation becomes the zero-temperature 
wave equation of a massive scalar in AdSs. The retarded solution to this equation is 

where Hu^^ is the Hankel function, d=A is the space-time dimension of the boundary theory, 
A/'q is an overall normalization, and 

^ = A-^. (5.9) 

The solution behaves like 

$ ~ Afgz'^-'^ (5.10) 

in the boundary limit J-O, where z = 2^fu. The divergence of f l5.10p as 2 — )■ for A > 
(i.e. vn} > 0) reflects the renormalization required of the corresponding operators in the 
4-dimensional gauge theory. In our discussion, we will be able to ignore the details of 
holographic renormalization prescriptions and simply summarize that (15. 8 p should approach 
z'^'^ (ph{q) as z approaches the (regulated) boundary, where 0b is the (renormalized) bound- 
ary source. Choosing (pb{q) = 1 defines the bulk-to-boundary propagator, which corresponds 
to (EHD with Afg ~ 1.1° 

The Hankel function goes through many oscillations, and so is well approximated by the 
geometric optics limit, when its argument is large compared to both 1 and u. In our case, 
this condition is parametrically equivalent to 

max(l, {RmY) 

U > Wmatch = 5 , (5.11) 



which may be also be written as 



M > max( ^— , Mmin). (5.12) 



-g2 



This generalizes the condition u ^ l/(— g^) previously discussed for the massless case. If 
we convolve (15. 8p with a high-energy source (11. 3p . we will not be able to use the particle 



In more detail, follow Ref. [13] and normalize the bulk-to-boundary propagator to be Zg ^ at z=z-b, 
where zb is arbitrarily small. Then 



Af, = z^ 



_ -1 

-A ' " 



ITT 



(iv^)^(4uB)'^/4ijW(y^iW) 



1 + 0{zb), 



and one takes zb — ^ at the very end of the calculation. In yet more detail, a systematic method for 



holographic renormalization is described in Refs. Ill4l3|. 
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approximation to figure out the details of wliat is liappening at m ~ Wmin, but we will be 
able to use it when the resulting wave packet propagates to m ~ provided m* ^ Wmatch, 
and so we may then use the particle approximation to calculate the stopping distance. 
Because we get the same stopping distance as for the massless case, we can take over 

.3 c f\q?\"\,c (E^'" 



The geometric optics approximation at m ~ m^, (and so this result for the stopping distance) 
will fail unless Wmatch ^ u^,. Using (13. ip and (15.111) . that condition requires 

^ Uax(l, (/?m)^) J • ^'-^'^ 

One might suspect that the particle approximation breaks down at the maximum possible 
stopping distance, in which case (I5.13P then gives that maximum to be 



X 



max{l, Rm) J \A 



1 . (5.15) 



We give a more detailed argument for this result in Appendix [Dl 

Eq. (15.151) implies that the maximum stopping distance decreases as the conformal di- 
mension A of the BPS source operator is increased. This qualitative feature is not novel to 
the strongly-coupled theory: it is true for the weaA;/?/-coupled theory as well. For the BPS 
operators, large A corresponds to an operator with roughly A powers of scalar fields, such 
as tr((/)^), where is one of the three complex scalar fields in the theory. In weak coupling, 
if we inject total energy E with such an operator, it will typically generate A particles that 
each have energy of order E/A. In weak coupling, the stopping distance of a particle with 
energy E scales as E^^'^ (up to logarithms), and so the stopping distance of the A particles 
each with energy E/A will scale as x^^^, ~ [E/ Ay/'^. 



B. Generalizing the power-law tail 

In this section, we investigate how the (x^)~^ power-law tail in Fig. [T] generalizes to other 
choices of source operator. We will take the source term in the gauge-theory action to be of 
the form 

e*^-^'AL(x)0(x), (5.16) 

where 0{x) is a scalar BPS operator with dimension A. We will see, however, that our 
result also applies to the case of O being a transverse-polarized R current, which was the 
case discussed in section [IVl 

Since the dependence (15.131) of the stopping distance on is the same as in the massless 
case, the only significant qualitative difference in the distribution of stopping distances will 
come from the distribution V{q^) of created by the source operator. The shape of this 
distribution is determined by the dimension A of the source operator, as we now describe. As 
discussed in Refs. j2|, 0]) temperature does not affect the initial creation of the jet, and so we 
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can simplify the analysis by evaluating V{q^) at zero temperature. Consider the probability 
density in q associated with a source operator O acting on the vacuum: 



5^(any|[0(g)]t|vac)*(any|[0(g')]Vac) 

= (vac|0(g) [0(g')]Vac) = iGy{q) {27rY5^'^\q - g'), (5.17) 



any 



where G> is the Wightman correlator of O. At zero temperature, it is related to the spectral 
density p of the operator by^^ 

zG>(g) = %°)p(g). (5.18) 

The distribution of jet 4-momenta is therefore given by the spectral density p{q)- At zero 
temperature, Lorentz and scaling invariance allow us to use simple dimensional analysis to 
know how p scales with g:^^ 

P(g) oc e(-g^) (-g^, (5.19) 

where = A — l(i as in (15. 9p . The 0{—q^) appears because only sources with time-like g^ 
produce persistent excitations at zero temperature. 

If O were a vector operator V^, like an R current, the dimensional analysis would be 
complicated by the fact that one could get factors of g^ associated with the vector index 
(rather than only factors of the virtuality g^). However the transverse spatial components 
q"*" of q^ vanish, by definition. So this complication does not arise for the transverse R 
current operator j-*- that we discussed earlier, and f l5.19p can also be used in that case. 

So far, we have only looked at the operator O and not the other factors in the source 
term (15.161) . We can rewrite (15.161) in g space as 

[~Al{q-k)0{q). (5.20) 

Jq 

Correspondingly attaching a factor of A2(g — k) to each 0{q) in (I5.17p . we get a probability 
distribution for g proportional to 

Oiq') p{q) \AL{q -k)? ^ ^(V) (-g')IA^(g - k)\\ (5.21) 

For a transverse-translation invariant source, the relative probability distribution for creating 
a jet with a given g+, where g^ ~ 4£^g+, is then 







J rfAg_ 


A 


f(?-+,Ag_)| 


2 


J dq+ dAq_ e{-q+) 




1 ^ 1 


Af(g+,Ag_)|' 



which generalizes (14. 7p . We will repackage this as 

2(3+v)/2 r 

n(g+) = -J^TT^ ^(-?+) k+^r ^.(9+^), (5.23) 



Alternatively, we could use the finite-temperature relation iGy{q) = [1 +n{q^)]p[q), where niuj) = 
[e^'^ — is the Bose distribution, and then use the fact that ~ £' ^ T in our problem. 

To get the dimension of p((?), use (|5.17p and note that 0{x) having dimension A means that the Fourier 
transform 0{q) has dimension A — d. 
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where 



T{^)jdq. |A?)(g+,g_)P 



(5.24) 



has been normahzed so that \I/;^(0) = 1 in the case of a Gaussian source envelope (jl.Sp . 
Following (14. lip , the distribution of stopping distances is then 



So the power-law tail in Fig. [T]has generalized to (x^)"'-^"'"^^-' = (x^)^"^'^, as shown in Fig. HI 
For the transverse R current operator, A = 3, which recovers our previous result of {x^)~^ 
in that case. 

C. Gaussian Source Envelope 

Throughout this paper, we have discussed two different scales x^^^-^^^-^ ~ {ELf/^ and 
^max ~ E^^^ characterizing the stopping distance. In f lS.lSp . we generalized the latter to 
^max ~ {E/ ^y/^ for the case of large A. Now we will discuss the similar generalization of 
^typical- -'■^ general, we will still have 



but the relation between the typical g+ of jets and the source envelope size L for large A will 
depend on details of how the source envelope Ai(g) falls off for large g+. That's because the 
probability distribution fl5.22p for the g+ of the jet involves a competition between the l^+l*^ 
factor which favors large |g+| and the |A2,(g+, Ag„)p factor which suppresses 3> . 
The typical value of g+ represents a balance between the two and will scale with v. For the 
sake of a concrete example, we consider here the case of a Gaussian source envelope (11.51) . 
In this case, the function \I'^(g+L) in fl5.23p and fl5.25p is simply 




(5.25) 




(5.26) 



(5.27) 




(5.29) 



(5.28) 
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Alternatively, one could compute the average value of x% 



stop 



directly from ()5.25p . giving x\. 



avg 
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VI. CONCLUSION 



The 5-dimensional particle picture provides a relatively easy way of understanding (from 
the gravity side of the calculation) the appearance of the scale {ELY^^ in jets created by 
finite-size sources in strongly coupled A/'=4 super Yang Mills plasmas. By thinking about 
sources with different types of momentum distributions, such as Figs. and b, we have 
learned that it is the range of of the source which determines the range of stopping 
distances. Making larger causes the jet to stop sooner. For any finite size L of source, 
the uncertainty principle implies that there will be a spread in the components of q of at 
least order 1/L and so a spread in of at least order E/L. As a result, almost all of the 
jets produced will travel distances < (ELY^^ [of order (ELY^^ in the case of Fig. [2^ and 
< {ELY^^ in the case of Fig. |2]d]. Events where a jet travels further (up to E^^^) will always 
be rare if the source size L is small compared to the maximum stopping distance scale E^^^. 

This interpretation, based on the 5-dimensional particle picture, provides an important 
clarification to our original derivation of the average distribution of charge deposition shown 
in Fig. [TJ This average includes an average over all events. From the original result, it was 
unclear whether on not Fig.[T]qualitatively tracks how the jet deposits its energy, momentum, 
and so forth on an event-by- event basis. It might have been that every single jet produced 
deposits some of its energy at ~ E^^^ and most of its energy at ~ {ELY^^- The 
success of the 5-dimensional particle interpretation, and in particular the success of (11. 8p . 
indicates that Fig. [1] instead reflects a probability distribution for how far the jet travels, 
and each individual jet dumps its energy and charge in a very localized region of x^. It 
should be possible to independently verify this conclusion by calculating correlations of the 
charge deposition at different distances, which we will leave to future work. 

Our analysis of massive 5-dimensional fields indicated that the maximum stopping dis- 
tance E^/^ depends on the type of operator used to create the jet and that the distance 
decreases as the conformal dimensions of that operator increases. That is, the maximum 
stopping distance depends on the type of high-energy excitation created. This may shed 
some light on a discrepancy between (i) the stopping distances found here and in Refs. 
j2|, 0|, which find a^max ~ E^^^, and (ii) those based on the evolution of classical strings in 
5 dimensions [l|, Isj, which find the parametrically smaller result x^^^ ~ {E / y/\Y^'^ , where 
A = NcQ^ is the large 't Hooft coupling. The gauge theory states corresponding to classical 
strings in the gravity dual may simply be states that are much more easily stopped by the 
strongly-coupled quark-gluon plasma than are the states created by the source operators 
considered in this paper. One could then ponder which (if either) might be more instructive 
for lessons about the theory of real QCD plasmas. In order to further clarify the differences 
between the two approaches, it would be interesting to find a 4-dimensional gauge-theory 
description of a source that could be precisely linked through duality to the 5-dimensional 
initial classical string configurations that have been used to study jet quenching. 
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Appendix A: Stopping distances in different dimensions 



The only relevant difference between AdSs-Scliwarzschild space and AdS^+i- 
Schwarzschild space is that the blackening function / = 1 — in the metric (I2.2p is replaced 
by / = 1 — v!^!"^ [14]. The stopping distance f l2.3p then becomes 



•^stop 



du 1'°° du 



Cd (WW [E\ ^^^^ 



where q = \/2 3(1+^, |~^)) B is the Beta function, and we have assumed d > 2. Taking e ~ 
l/L, the dominant stopping distance (II .Qp for our original source of Fig. then generalizes 
from {ELy/^ to 

Epical ~ {ELY'-'^/''. (A2) 

The quickest way to estimate the maximum stopping distance, generalizing E^^^, is to esti- 
mate when the stopping distance ( 1A2P becomes as large as the source itself, as we did for 
d=4: in section [11 The result is 



.3 

^max 



The last result can also be obtained from a wave analysis by analyzing the poles of the 
retarded bulk-to-boundary propagator, just as was done for d=4: in Ref. The scale of 
the exponential decay in fll.lbp was determined by the imaginary part of the propagator 



pole closest to the real axis. Here we simply follow section 4.6.2 of Ref. generalizing to 
arbitrary d. The massless field equation for A± is 



,2 4:Eq+ - u'^/^E^ 



u 



= 0. (A4) 



Changing variables to ^7 = e *27r/(d+2)^4/(c(+2)^^ ^^le field equation becomes 

'-dfj+{u^'-'^/'-^)]A^ = (A5) 

where 

Poles of the bulk-to-boundary propagator occur when the Schrodinger-like equation (lASp 
has a zero-energy bound state that vanishes at the origin. The smallest value of a for which 
this occurs is 0(1), from which flA6P gives that the pole closest to the origin has 

Imgf^~E-('^-2)/('^+2). (A7) 

So the response to the source falls exponentially as 

for large x"*", for some constant k^. This behavior is consistent with flA3p . 
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Appendix B: Null geodesic in AdSs-Schwarzschild 

For the sake of keeping this paper self-contained, we give here a brief derivation of (12.11) . 
A null geodesic has 

= {dsY = dx^ dx" + dx^ g^s dx^, (Bl) 

and so 



dx^ 1 



d\ Jg55 



dx^ dx^ 



1/2 



d\ d\ 



(B2) 



where we will take A to be any affine parameter for the trajectory. Because of 4-dimensional 
translation invariance, 

(B3) 

is conserved and proportional to g^, so that 

dT^ 

oc g^y.. (B4) 



d\ 

Dividing this equation by (1B2I) gives 



dx^" ^ — gf'^q^ 

955 I „ M/2 ' 



dx^ ^ '''''' {-q.g'^^qpY/ 
which in turn gives (12. ip . 

Appendix C: More on the geometric optics approximation 

Here, we will go into a little more detail about the conditions for the geometric optics 
approximation, for space-time backgrounds with 4-dimensional Poincare invariance. For 
the sake of concreteness, we will consider the case of the source (14. ip used in our earlier 
work, where the source operator is a transverse-polarized R current and so is dual to a 
5-dimensional transverse vector field A±_ = e^A^ in the gravity description. The equation of 
motion for v4x(g, x^), where q is the 4-momentum, has the form 

^ : d,{^ g^^g^'d.A^) = g^^q.g^'q^A^ (CI) 



where g'^^ is the component of the inverse metric in the direction of the polarization, e.g. 
g-^^ = g^^ = g"^"^ for q in the x^ direction. Now switch to coordinate 

(C2) 



which parametrizes proper length in the direction of the fifth dimension, and note for future 
reference that 

de = ^ds. (C3) 



The equation of motion is then 

w-'de{wdeA^) = q.g^'q.A^. (C4) 
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where 

w = ^-g{A) 9^^ (C5) 

and (7(4) = g^^g is the determinant of the 4-dimensional part of the 5-dimensional metric. 
(As another simple example, one could consider the case of a source operator dual to a 
massless scalar field, which would corresponds instead to taking w = ^—(7(4) in the equation 
of motion.) Now define 

a = y/wAi_ (C6) 



to get 



(C7) 



This looks just like a one-dimensional quantum mechanics problem in i with wavenumber 
k{i) = ^2m{E - V{i)) replaced by 



The WKB condition that the wavelength X{i) = 271 /k{i) in such a quantum mechanics 
problem not change significantly over distances of one wavelength is deX{i) ^ 1. 
For the metric (Q, w = {R/2fp/^u-^, and (JCH]) becomes 



k{u) = - Q') + (C9) 

with the help of ( 1C3|) . For the m's of interest to our discussion of particles, which are 
M ~ ^ 1, the (1 — 2u^)/4:f term under the square root is negligible, giving 

Ku) ^ ^^qilg^u ^ '^^u{uW-q')- (CIO) 

Again using (IC3I) . the condition di\{t) = de{2TT/k) <^ 1 then gives (13. 6 p for u ~ m*. 



Appendix D: Maximum Stopping Distance for High-Dimension Source Operators 

In section IV At we saw that the particle picture breaks down when > x^^^ with 

3 f E V'"" [eV'"' , , 

^-^ Uax(l,i?m) J ~UJ ' ^ 

and we suggested that this x\^ was the furthest that jets would propagate — that is, that 
energy or charge deposition at larger distances would be exponentially suppressed, similar 
to our previous A=3 result of (ll.lbl) . In this appendix, we will discuss how (1D1I) arises in a 
wave analysis following the methods of Ref. We will focus on the case of A ^ 1. As in 
the main text, we assume that the large E limit is taken jirst^ and only then do we consider 
large A. 

There is a subtlety to the results we will find. In section ID 2\ we will analyze the expo- 
nential fall-off of jet charge deposition at very large x by finding the location in the complex 
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g+ plane corresponding to the first quasi- normal mode of the bulk field. For large A, we will 
find an exponential fall-off of the form 

/ 2c' \ 

-p(-:^T77Ai7^j> m 

where c'^ is a constant. This is the large-A version of the exponential in (ll.lbl) . The scale 
that determines the rate of exponential fall-off in (lD2p is £'^/^/A^/^. One might naturally 
guess that exponential suppression therefore applies whenever x^ ^ E^^^ / A^^^ and so guess 
that x^jjx ~ ii^^/^/A^/^ instead of fIDll) . This guess fails, however, for reasons we shall now 
outline. 



1. Overview 



To understand the issues involved, we briefly highlight some relevant aspects of the m=0 
calculation from Ref. The main part of the calculation there involved computing the 
bulk response A to the high-energy source on the boundary, given by 

A{x, u)= I g\q, u) ~AL{q - k) e*^■^ (D3) 

Jq 

where Q is the bulk-to-boundary propagator. Our result for the stopping distance came 
from extracting the behavior of A near the horizon, u — )■ 1.^'^ The critical part of the q 
integration was the integral over For x^ <^ E^^^ (in the m=0 case), we found that 
we could deform the g+ integration contour in the complex plane so that the integral was 
everywhere exponentially suppressed except at a saddle point 

of (13. 2p . This contour is depicted in Fig. El Parametrically far into the interior of the shaded 
region indicates places where the magnitude of the integrand is exponentially suppressed. 
Parametrically far into the unshaded regions indicates places where it is exponentially large. 
The dashed line depicts a line of poles of the bulk-to-boundary propagator corresponding 
to quasi-normal modes. In the WKB approximation to that propagator, this line of poles 
became a cut. 



Readers may wonder at the connection between (i) studying u — )■ 1 and (ii) studying m ~ ^ 1 as 
in the particle arguments in the main text of this paper. The point is that how far the particle travels 
is determined by where it is when u — > 1, but the integral which gives that distance is dominated by 
u ^ U-): <^ \. In the wave analysis of Ref. [4], we studied the response at late times, corresponding to 
u— >1 for the bulk excitation. Our results for the near-horizon bulk response were determined by the 
term in the WKB exponent S for Q (see eq. (4.51) of Ref. j^). But this term was generated 
by the region of the integral that gave S. See, for example, eq. (D13) of Ref. 3, which is proportional 

to the current paper's particle stopping distance integral (12. 3p . 

In Ref. we expressed formulas in terms of X+ = a;+ — t(m) instead of x^. As discussed in that paper, 
the late-time response is localized to x~ ~ ^ '''(w) (see eq. (4.49) of Ref. 1^]), and so ~ 2x^. 
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FIG. 6: Integration contour in the g+ complex plane for saddle point approximations to the g+ 
integral in ()D3p . This is a slightly simplified version, appropriate for n— t-I, of Fig. 13b of Ref. 
0]. The location of the saddle point is marked by the large dot. The circles indicate different 
parametric scales for \q+\. 




For ^ E^/^, the regions of exponential suppression for the integrand are shown in 
Fig. [7^. Saddle point methods are unreliable. Instead, choose the integration contour 
shown there. The piece that goes around the line of poles can be re-expressed as a sum 
of contributions from each pole, as shown in Fig. [7)d, which depicts a magnification of the 
neighborhood of the origin of Fig. [7^. The e*''+'^^ piece of the e*"^'^ in (1D3P causes the 
contributions from the poles to be exponentially suppressed according to their distance 
Img+ from the real axis. For ^ E^^^, the nearest pole dominates and produces the 
exponential fall-off ( ]l.lb|) of the jet's charge deposition. 
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FIG. 8: Like Fig. [6] but for the case of large A and ^ E^^^ / A^^^ . The cross-hatched region 
represents < A'^/^/E^/^. 



Now we return to the massive case, with A ^ 1, and ask what happens for 



A4/3 



< a;^ < 



AV3- 



(D5) 



This is the interesting case, where (i) the geometric optics and particle arguments of section 
IVl indicate that the charge deposition is not exponentially suppressed but (ii) the guess we 
might make based on flD2p suggested that it is suppressed. As we shall discuss in section 
ID 31 the mass does not significantly affect the massless picture of Fig. [6] provided 



A4/3 



(D6) 



which is the condition fl5.14p discussed in the main text. For <C E^^^/A^^^ as in fIDSp . 
the condition flD6P is satisfied at the saddle point flD4p and for larger So, as long as 
we are careful to route the contour as in Fig. [HI we can take over the methods of the m=0 
calculation, make a saddle point approximation to the g+ integral (which corresponds to 
making the geometric optics approximation), and so find a result that is not exponentially 
suppressed. 

But now consider a large calculation, along the lines of Fig. [71 We shall see in section 
ID 21 that the closest pole to the real axis has 



Im q: 



pole 



A4/3 



(D7) 



That is, the line of poles in Fig. [H] begins at the edge of the hatched circle. For x^ in the 
range of ( 1D5I1 . picking up the poles is not so useful. Fig. [9] shows the large- A analog of Fig. [7] 
for this x^ range. As one looks at poles progressively further from the origin, the exp{iq ■ x) 
factors leads to suppression, as before, but the factor grows exponentially (as well as 
oscillates), and this exponential growth dominates the integrand in the unshaded region of 
Fig. [9] So one cannot approximate the integral by the contribution from the nearest pole, 
and asymptotic formulas like (ID2p do not apply in this case. 
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Since the saddle point approximation is controlled and gives a result that is not expo- 
nentially suppressed in the range fIDSp . whereas the approximations that lead to flD2p are 
not valid there, we conclude that the maximum stopping distance is of order E^^^ / A^^^ and 
not EV3/a4/3_ 

What happens as one continues to increase For ^ E'^/^/A^/^, the saddle point 
approximation breaks down and Fig. [9] becomes Fig. [TOl In this case, the calculation is 
dominated by the first pole and (lD2p applies. We will not attempt here to calculate the 
details of the transitional behavior at a; ~ xl,„^}^ 



However, for the sake of showing that (jD2p is mathematically consistent with x^^^^ ~ E^^^ / A^^"^ , we 
point out that {exp[2^jy|— ^^5^] + 1}"^ is an example of a function that is unsuppressed for x ^ a^max 
but decays like (jD2p for x x^^^^. 
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2. Quasi-normal modes 



To find the poles of in the complex g+ plane for large E, we follow the method of Ref. 
For the sake of simplicity, we will focus on the case of a massive bulk scalar field, whose 
equation of motion is f l5.7p . As discussed in Ref. the pole positions at high energy are 
determined (up to small corrections) by the nature of the equation of motion for u <^ 1. It 
will be more convenient to work with the variable z = 2y/u instead of u. Writing 



z 



(d-l)/2 



the equation of motion for u ^ 1 becomes the Schrodinger-like equation 
with potential 



[Rmf + 



d?-i 



Taking the high energy limit and and setting (i=4, this is 



with 



and 



Viz) 



15 
4 



= {Rmf + 

Following Ref. make the change of variables from u to 



(D8) 
(D9) 
(DIO) 

(Dll) 
(D12) 
(D13) 

(DM) 



which turns the retarded boundary condition at large u into the requirement that ip be real 
and exponentially falling. In terms of z^ this redefinition is 



./6^1/3. 



The resulting equation is^^ 
with 



V{Z) ^ - 



Z2 



(D15) 
(D16) 

(D17) 



For comparison with Ref. [3], one may write a similar equation in terms of U by defining (p ~ u''' 
giving [-9?r + U 
{Rmf 



iE^nilJLl - i],/, = for d = 4. This reduces to (4.66) of Ref. [41 for A±_ in the case 
3. The formulation in the current paper in terms of Z is more 



c/ ' ^ ' HP u 
3, corresponding to A 

convenient because the pole locations q+ can be identified as proportional to the bound state energies of 
a Schrodinger potential V{Z). Note also that M.^ plays a roll analogous to angular momentum squared in 
the Schrodinger problem (jD16l) . with jZ"^ like a centrifugal potential and the large- limit analogous 
to a large angular momentum limit. 
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and a defined in terms of g+ as in Ref. 

Solving ( 1D16P witli tlie desired boundary conditions is equivalent to setting a to be twice 
the bound-state energies associated with the potential V{z). For M. ^ 1, these can be well 
approximated by treating V{z) in harmonic-oscillator approximation around its minimum. 
The result is 

for n = 0, 1, 2, ■ ■ ■ }^ Using fIDlSp . we find that the first pole in g+ is a distance of order 
/\4/3^^i/3 ^Yom the real axis, but the spacing between successive poles in Fig. [TOb is only of 
order A^^^ / E^^^ . The specific result for oq determines 

_ VSoo ^ 3^3 
^i-4AV3- 16 ^^^^> 

for the exponential fall-off f lD2p in the case A ^ 1. 



3. Saddle point analysis 



In section ID 11 we claimed that a large mass m for the bulk field does not qualitatively 
change the massless saddle-point picture of Fig. [6] except inside the (avoidable) hatched 
region of Fig. [HI Here, we will briefly sketch why. For Rm ^ 1, the condition for the 
validity of the WKB approximation^^ is satisfied in the small u regime u <^ Umm as well 
as in the oscillatory regime u ^ Umin- (The turning point u ~ u^ain can be avoided simply 
by analytically continuing around it, as in the textbook discussion of WKB in Ref. [HI). So 
we may use WKB all the way to the boundary u = ub- 



oc e'^ = exp 



du' qsiu') 



(D21) 



where, for the sake of simplicity of presentation, we will suppress showing the WKB prefac- 
tor. For the massive case, 



1 /u^lgrp — g2 (RrnVf 

^^(") = 7V — u 



The integrand in flD3p then has exponential dependence 

and the saddle point of its integral is determined by 

dS d 







M L 



q ■ X + dx qs{x 



(D22) 



(D23) 



(D24) 



Here we label the first pole qq. In Ref. [4] we instead called it oi. 

Essentially: that the derivative of the WKB exponent does not change significantly over one e-folding or 
oscillation. 
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which gives 

= - f dx^^. (D25) 
J dq^ 

Together with (1D22I) . this simply reproduces the particle-based formula (15. 3p for the stopping 
distance. Approximating c::^ E and solving for g+ in terms of x^ will then give the saddle 
point q^ for the g+ integration in fIDSp . We've already discussed the effect of the mass on the 
particle stopping formula (15. 3 p back in section IV A[ Tracing the discussion of section IV Al 
backward, x^ <^ {E/ Ay^^ corresponds to —q+ given by e 2> {E / A^)~^^^ , which corresponds 
in turn to Umin ^ u^,. That's precisely the case where the mass had a negligible effect on 
the relationship between x^ and g+. In consequence, the mass m will not have a significant 
effect on the determination of the saddle point q*^_ for x^ ^ (E/ Ay^^. 

What about the behavior of the integrand elsewhere along the contour in Fig. [HP A 
discussion of the WKB exponent S of (]D2ip is complicated by the divergence (I5.10p of the 
bulk-to-boundary propagator on the boundary, which shows up as a logarithmic divergence 
[oc In Mb) of the integral in (lD2ip . We will briefly indicate in section [D 41 how one can do 



a WKB analysis that avoids this divergence, but such details lose the forest for the trees. 
More simply, the Iumb divergence of / du gs is independent of g+ and so does not affect 
the g+ dependence of the integrand in (lD3p . and so it will only affect the result by overall 
factors. To focus on the question of whether the mass makes a significant effect on the 
dependence, look at the effect of the mass on dS/dq^ instead of on S. So look at 

The first term is mass independent, and the second term is just once again our integral 
for the particle stopping distance as a function of g+, given by the right-hand side of (15. Sp . 
though with an imaginary part even for real negative g+ due to integrating over u < Mmin- 
However, if |g_|_| ^ A^^^ / E^^^, then the effect of the mass on this integral will have negligible 
relative magnitude, just as in the previous discussion concerning the location of the saddle 
point. 



4. Avoiding WKB exponent divergences 



Finally, we sketch how one could set up a finite WKB integral if one wanted to carry 
through the analysis of this appendix in more detail than we have given. To get the nor- 
malization (I5.10p appropriate for the bulk-to-boundary propagator Q, we want $ to give 



u 



{d-A)/2 
B 



at m=mb. So 



(d-A)/2 

g exp 



du'q^{u') 



(D27) 



where we have again suppressed showing the WKB prefactor, other than the overall power 
of Mb. We can trade the divergent ^^^'^ normalization factor for a finite m('^~'^)/2 by 
introducing a compensating change in the exponent: 



7R 



u 



{d-A)/2 



exp 



du' 



q5[U 



2u' 



(D28) 
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In the large A limit (required for our WKB analysis in the region u <^ Wmin), A — c? ~ Rm, 
and so we will replace the last equation by^° 



u 



I /(w 



1 /M'2|q|2_g2 (Rm)^f{u') Mm 



u' 



/2 



2u' 



(D29) 



Now the integral in the exponent is finite if we take the limit u-q — )■ 0, and so the ap- 
propriate WKB expression (still suppressing showing the original WKB prefactor) can be 
approximated as 



gE^„(<i-A)/2j,^p 



1 /u'2|g|2-g2 {RmYf{u') .Rm 



u' 



/2 



2u' 



(D30) 



One may then use this WKB formula to pursue a more detailed analysis. The —iRm/2u' 
subtraction in the integral cancels the original integrand for u -C Umin and so keeps the 
integral finite. Its contribution for u ^ Umm will introduce a additive piece of approximately 
|i?mlnMmin — |AlnMmin iu the WKB exponent S. That corresponds to a multiplicative 

factor of WjnirL^^ i'^ result for e**^ and so -u"^ in \A\'^. Using (15. 5p for and then the 
saddle-point value ( ]D4p for g+, this factor is 



^, (D31) 
which is just the dependence of the power-law tail on z/ ~ A that we previously found in 

(ESS}. 
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